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SUMMARY 

Probabilistic assessment of the maximum blade response of a mistuned 
rotor (bladed disk assembly) is performed using the computer code NESSUS. The 
uncertainties in natural frequency, excitation frequency, amplitude of excita- 
tion and damping have been included to obtain the cumulative distribution func 
tion (CDF) of blade responses. Advanced mean value first order analysis is 
used to compute the CDF. The sensitivities of different random variables are 
Identified. The effect of the number of blades on the mistunuing of rotor is 
evaluated. It is shown that the uncertainties associated with the forcing 
function parameters have significant effect on the response distribution of 
the bladed rotor. 


INTRODUCTION 

This paper demonstrates the applicability of the computer code NESSUS 
(Numerical Evaluation of Stochastic Structures Under Stress) developed at NASA 
Lewis Research Center (refs. 1 and 2) to obtain probabilistic assessment of a 
mistuned bladed disk assembly. The inherent uncertainties in the blade geome- 
try, material properties, damping, mass, etc. result in uncertainties of 
dynamic characteristics. Due to these uncertainties, each blade on a rotor 
disk does not have same fundamental frequencies and damping characteri sti cs . 
These differences are referred to as mistuning. When the blades with such 
small variations in frequencies interact together on a rotor disk, their indi- 
vidual amplitude and dynamic stresses are altered significantly. This is noth 
ing but the effect of mistuning. The overall disk assembly will have many 
vibration modes which are very close to each other. The closely spaced modes 
amplifies the responses significantly. Many researchers have investigated the 
effect of mistuning both in a deterministic (ref. 3) and nondetermi ni sti c (pro 
babilistic) (ref. 4) sense. Main focus of those investigations has been on 
uncertainties in the blade characteristics. However, in reality the uncertain 
ties are also associated with the loads. 

This paper Includes both structural and load uncertainties. The nondeter 
ministic analysis can be performed by Monte Carlo simulation technique or any 
other advanced probability integrator method. Monte Carlo simulation 



technique is the most reliable method but it is computationally inefficient 
and expensive. The probabilistic methods available in NESSUS have been 
applied to the problem of mistunlng and found to be efficient and economical 
to compute even the low level probability responses. The unique feature of 
NESSUS is that the reliability analysis can be performed directly for a given 
statistical distribution of resistance. The outcome of the reliability analy- 
sis can be easily and directly used for the design and risk analysis purposes. 


NESSUS - BRIEF DESCRIPTION 

NESSUS is a comprehensive Probabilistic Finite Element Analysis computer 
code to perform static, dynamic, buckling and nonlinear analysis. Various pro- 
babilistic and finite element methods are integrated into NESSUS in a modular 
fashion to solve any general purpose structural analysis problem. Different 
modules in NESSUS can be used independent of each other. Also, all the 
modules have been linked to automate the reliability analysis of any structure 
from beginning to end. Three basic major modules namely; NESSUS/PRE, 

NESSUS/FEM and/NESSUS/FPI form the core of code. The usage of NESSUS in 
general involves the identification of primitive variables, their statistical 
distributions, structural geometry, loads, boundary conditions, etc. Any 
strength of statistical correlation can be input to NESSUS. 

NESSUS/PRE is a preprocessor which allows the user to describe any 
Gaussian random field over the spatial domain defined in terms of discrete 
finite element nodes. The Gaussian random field is defined in terms of random 
variables at nodes with their mean, standard deviation and appropriate form of 
correlation. Correlated random variables are decomposed Into a set of uncorre- 
lated Independent vectors using modal analysis by NESSUS/PRE. 

NESSUS/FEM is a general purpose finite element code to perform structural 
analysis and evaluation of sensitivity due to variation in different uncorre- 
lated primitive variables. The response surface defined in terms of random 
variables required for probabilistic analysis in NESSUS/FPI is obtained dis- 
cretely in NESSUS/FEM by perturbing Independent random variables obtained from 
NESSUS/PRE. NESSUS/FEM incorporates an efficient perturbation algorithm to 
compute the sensitivity of random variables. The response sensitivity data 
generated is stored on a file using linked list database structure. NESSUS/FPI 
has several advanced reliability methods including fast Monte-Carlo simulation, 
the most efficient one being the Fast Probability Integrator. This module 
extracts the data for a given response from the database created by NESSUS/FEM. 
Since FPI needs explicit response function, the extracted numerical data is 
used to develop performance model in terms of Independent random variables. 
Using thus developed performance model and statistical distributions of primi- 
tive variables, the cumulative distribution is obtained for a given response 
using fast probability Integration. 


BLADED DISK ASSEMBLY PROBLEM 

The objective of this paper is to demonstrate the applicability of NESSUS 
to the problem of a bladed disk assembly as described in the introduction. 

The aim is: (1) to obtain cumulative probability distribution function (CDF) 
of the blade amplitudes and (2) to evaluate the effects of: (a) number of 
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blades on a disk (b) the statistical distribution of primitive variables and 
(c) the order of excitation. The following cases are simulated to achieve 
these objectives: 

(1) Disk with 10 blades 

(a) First order excitation 

(b) Second order excitation 

(c) Mean excitation frequency smaller than mean blade frequency 

(2) D1 sk wi th 20 blades 

(a) First order excitation 

(b) Second order excitation 


(A) Finite Element Model 

The finite element model consists of geometry and loads. 

Geometry . - A lumped mass-spring model representing the blades and rotor 
disk has been used for the purpose of analysis. The disk is represented by 
lumped masses connected together by a set of radial and circumferential 
springs. The masses representing the blades are mounted on each disk lumped 
mass as shown in figure 1. Each blade is modeled by two lumped masses one of 
which corresponds to the aerofoil region and the other to the shank and 
fur-tree region. The lumped masses are connected by beam elements. The 
radial and circumferential beam elements for the disk have only extensional 
stiffness representing the respective stiffnesses for the disk. The beam ele- 
ments connecting blade masses have only bending stiffness to represent the 
flexural mode of vibration of the blade. Therefore, it is important to note 
that only the bending blade mode of vibration is considered in the analysis 
since It normally dominates the responses. Thus, the dynamic interaction 
between blades comes through disk springs and masses. 

Loads . - The force induced on blades is due to variations in the flow 
around circumferential direction. The variations in the flow are considered 
to be periodic in nature. Therefore, the excitation would be harmonic which 
is characterized by amplitude, frequency of excitation and phase angle. The 
magnitude of amplitude does not matter in case of linear analysis, however its 
statistical distribution does since the assessment is made on the responses 
normalized to those of tuned rotor. The frequency of excitation is selected 
as Integer multiple of engine rpm. Since, the excitation is harmonic in 
nature, it differs from blade to blade in terms of phase angle only. As the 
blades are fixed around the circumference, the phase difference, 4» between 
blade's excitation is given by 


A 2Mir 

♦ = nr 

where M is the integer multiple between engine rpm and frequency of 
excitation 


N is the number of blades on a disk 
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The modal damping used in the analysis Is 3 percent (includes material, 
structural and aerodynamic). 


(B) PROBABILISTIC MODEL 

The following variables are considered to be random in the analysis: 

(1 ) Frequency of blade 

(2) Damping 

(3) Frequency of excitation 

(4) Amplitude of excitation 

The frequency of the individual blade depends on Its geometry, material 
property and the mass. Since frequency is not a primitive variable input to 
NESSUS, modulus of elasticity, E of blade is substituted for it as a random 
variable. This is done in such a way that the mean value of E gives mean 
frequency and one standard deviation (SD) of E from its mean gives frequency 
equal to one SD from its mean. In this simulation the uncertainties in the 
geometry and the mass are represented by equivalent uncertain moduli. 

Normally, the rotor disk is balanced by placing blades of the same moment 
weight In opposite slots. To Include such a physical requirement In a probabi- 
listic model, it is assumed that the blades in opposite slots are fully corre- 
lated statistically. The mean, standard deviation and distributions of 
different random variables are given in table I. 

Since, NESSUS/FPI requires a performance model (response surface) of 
response to compute the CDF, the sensitivity of each random variable is 
obtained by perturbing them one at a time. Each blade frequency (modulus of 
elasticity) Is considered to be an independent random variable in the analysis 
to develop a suitable performance model of mistuning. 


(C) MATHEMATICAL ALGORITHM DISCUSSION 
The governing equation of motion for the bladed disk is given by 

[M] {X} + [C] {X} + CK3 {X} = {Ae 1 <“*-♦>} 

where X, X, X is acceleration, velocity and displacement respectively 

M mass 
C damping 
K stiffness 

A amplitude of excitation 
b) forcing frequency 
4> phase angle 
t time 

NESSUS uses the modal response transfer function approach to calculate 
the modal responses. Modal responses are superposed to obtain total response. 
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One of the major problems encountered In the analysis of mistuned rotors 
Is to extract fairly large number of closely spaced eigenvalues. There are at 
least as many closely spaced frequencies as the number of blades. Also, the 
fact that the first blade on a disk Is connected to the N th blade, there 
exists a nonzero term In the corner of stiffness matrix. Thus the stiffness 
matrix becomes sparse rather than banded. NESSUS has a powerful profile stor- 
age scheme and a subspace Iteration eigenvalue extraction algorithm. When the 
number of closely spaced modes Increase, the tolerance provided for Iterative 
procedure should be small. The algorithm did not show any sign of numerical 
Instability and the eigenvalues were extracted In relatively few Iterations. 
The other salient feature of the algorithm. In conjunction with perturbation 
procedure. Is that the eigenvalues and eigenvectors of the previous perturba- 
tion are used as Initial trial vector for the current perturbation In order to 
obtain rapid convergence and Increased accuracy. 

Advanced mean value first order analysis (AMVFO) (ref. 5) has been used 
to obtain CDF of amplitudes. The procedure Involved can be summarized as 
fol lows: 

(1) Obtain the linear approximation of response, Z](5() based on perturba- 
tion about mean value. 

(2) Compute values of Independent variables, X* for selected CDF values 
(mean value first order analysis). 

(3) Recompute Z(X*> to correct the response for the same CDF In step (2) 
(Advanced mean value first order analysis). 

It has been found that the AMVFO gives very accurate results efficiently. 


DISCUSSION OF RESULTS 

Cumulative probability distribution of normalized amplitude of the maxi- 
mum blade response on a disk are plotted. The normalized amplitude Is defined 
as the ratio of maximum blade response of a mistuned rotor to that of the 
tuned rotor. The tuned rotor Is the one for which all the individual blade 
frequencies are the same. The normalized response indicates the amplification 
of response due to mlstunlng. Figures 2 and 3 depict the CDF of normalized 
blade amplitude for a rotor with 20 blades and figures 4, 5 and 6 for the ro- 
tor with 10 blades. As expected, the results show that for the same probabil- 
ity level for a given rotor, the amplification Is higher when the excitation 
frequency Is closer to the natural frequency of a blade. Also, it is 
observed that for a given order of excitation, the amplification increases as 
the number of blades Increase. This Is due to the fact that as the number of 
blades Increase, a higher number of vibration modes get closer to each other. 
The following Is very interesting and worth noting (refer to table II): 

(1) When mean excitation frequency is larger than the mean blade natural 
frequency, the uncertainties in the excitation frequency are more sensitive to 
response CDF than those of blade frequency. 
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(2) When mean excitation frequency Is smaller than the mean blade natural 
frequency, the uncertainties of blade frequency become significantly sensitive 
to response CDF as compared to those of excitation frequency. 

The sensitivity of the frequency of a blade or the frequency of excita- 
tion depends largely on the difference between the mean of each frequencies 
and their coefficient of variations. If: 

(1) The coefficient of variation for the blade frequency is small com- 
pared to that of the excitation frequency. 

(2) The means of the two frequencies are fairly far apart. 

(3) The mean excitation frequency Is larger than mean blade frequency. 

Then It is obvious that the uncertainty in the excitation frequency will 
dominate the response distribution of the rotor. The point to note Is that it 
Is very Important to consider uncertainties in the blade frequency and excita- 
tion forces simultaneously in order to properly assess their respective sensi- 
tivities on the response of the rotor. 


CONCLUDING REMARKS 

It Is demonstrated that the uncertainties of mlstuned bladed disks can be 
evaluated using NESSUS. Structural and modeling considerations related to the 
problem have been included. The effect of the number of blades, and the order 
of excitation have been Identified. The specific findings are as 
fol lows : 

(1) When the mean excitation frequency is larger than the mean blade 
frequency: 

The effects of uncertainties In the excitation frequency to mlstunlng is found 
to be very significant as compared to the effects of uncertainties In the 
blade frequency. Therefore, It Is very Important to simultaneously con- 
sider both structural and load uncertainties in the analysis. 

(2) When the mean excitation frequency Is smaller than the mean blade 
frequency: 

The effects of uncertainties In the blade frequency are more significant in 
comparison to those of excitation frequency. 

These findings, (1) and (2), and corresponding sensitivities help 
designers achieve robust rotor designs. For example, when the scatter in the 
blade frequency Is larger and uncontrollable, then the rotor should be designed 
to operate beyond resonance and vice versa. The blade can be designed for a 
desired probability level and Its corresponding amplification of response. 

Also, the sensitivities provide directions (through the primitive variables) 
to enhance quality control In order to reduce the uncertainty in the rotor 
response and consequently the risk. 
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TABLE I. 


[Statistical distribution of all primitaive variables considered 
i s normal . ] 


(a) Disk with 10 blades. 


Primitive variable 

Mean 

Coefficient of 
variati on , 
percent 

Frequency of blade (cps) 

959.0 

1.0 

Damping 

3 percent 

5.0 

(a) First order excitation 



Frequency (cps) 

1909.86 

5.0 

Ampl i tude 

100.0 

5.0 

(b) Second order excitation 



Frequency (cps) 

2864.8 

5.0 

Ampl i tude 

100.0 

5.0 

(c) Mean excitation frequency 
smaller than mean blade 
frequency 



Frequency (cps) 

477.5 

5.0 

Ampl i tude 

100.0 

5.0 


(b) Disk with 20 blades. 


Pri mi t i ve vari abl e 

Mean 

Coefficient of 
variati on , 
percent 

First of blade (cps) 

959.0 

1,0 

Damping 

3 percent 

5.0 

(a) First order excitation 



Frequency (cps) 

1909.86 

5.0 

Apl i tude 

100.0 

5.0 

(b) Second order excitation 



Frequency (cps) 

2864.8 

5.0 

Ampl i tude 

100.0 

5.0 


TABLE II. - SENSITIVITY FACTORS 
[Case (a), (b), (c) etc, refer to those of table I.] 


(a) Di sk wi th 10 bl ades . 


Probabi 1 i ty 

Primative Variable 

Sensi ti vi ty 

(a) 0.99999 

Blade frequency 

0.15 


Excitation frequency 

.916 

(b) .99999 

Blade frequency 

.017 


Excitation frequency 

.953 

(c) .99999 

Blade frequency 

.527 


Excitation frequency 

.657 


(b) Disk with 20 blades. 


Probabi 1 i ty 

Primative Variable 

Sensi ti vi ty 

(a) 0.99999 

Blade frequency 

0.124 


Excitation frequency 

.915 

(b) .99999 

Blade frequency 

.051 


Excitation frequency 

.91 
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CUMULATIVE PROBABILITY 


i TH BLADE 




FIGURE 1. - LUMPED MASS-SPRING MODEL FOR BLADED DISK. 


FIGURE 2. - MISTUNING OF BLADED DISKS. 20 BLADES 2nd 
ORDER EXCITATION. 



FIGURE 3. - MISTUNING OF BLADED DISKS. 20 BLADES 3rd 
ORDER EXCITATION. 


FIGURE 4. - MISTUNING OF BLADED DISKS. 10 BLADES 2nd 
ORDER EXCITATION. 



FIGURE 5. - MISTUNING OF BLADED DISKS. 10 BLADES 3rd 
ORDER EXCITATION. 
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